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Abstract 


Analyzing test particles falling into a Kerr black hole, we study gravitational 
waves in Brans-Dicke theory of gravity. First we consider a test particle plung¬ 
ing with a constant azimuthal angle into a rotating black hole and calculate 
the waveform and emitted energy of both scalar and tensor modes of grav¬ 
itational radiation. We find that the waveform as well as the energy of the 
scalar gravitational waves weakly depends on the rotation parameter of black 
hole a and on the azimuthal angle. 

Secondly, using a model of a non-spherical dust shell of test particles falling 
into a Kerr black hole, we study when the scalar modes dominate. When a 
black hole is rotating, the tensor modes do not vanish even for a “spherically 
symmetric” shell, instead a slightly oblate shell minimizes their energy but 
with non-zero finite value, which depends on Kerr parameter a. As a result, 
we find that the scalar modes dominate only for highly spherical collapse, 
but they never exceed the tensor modes unless the Brans-Dicke parameter 
‘.vbd ^ 750 for a/M = 0.99 or unless ujbd ^ 20,000 for a/M = 0.5, where M 
is mass of black hole. 

We conclude that the scalar gravitational waves with ujbd "S:, several thou¬ 
sands do not dominate except for very limited situations (observation from 
the face-on direction of a test particle falling into a Schwarzschild black hole 
or highly spherical dust shell collapse into a Kerr black hole). Therefore ob¬ 
servation of polarization is also required when we determine the theory of 
gravity by the observation of gravitational waves. 
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I. INTRODUCTION 


After the discovery of the binary pulsar PSR1916+13 by Hulse and Taylor |^, we be¬ 
lieve in the existence of gravitational waves emitted from compact objects. Their continual 
observations over 20 years show that both the observed change rate of the orbital period 
and that predicted by general relativity agree to within 0.3% accuracy p|. 

Currently, a number of worldwide projects for detecting gravitational radiation directly 
using a km-scale laser interferometer such as LIGO, VIRGO, GEO600, and TAMA are 
progressing. In particular, the LIGO project 0 will begin operating within a few years. 
Among the many possible sources of gravitational waves, one expected plausible source is a 
coalescing binary, composed of neutron star-neutron star (NS-NS), neutron star-black hole 
(NS-BH), or black hole-black hole (BH-BH). These direct detections of gravitational waves 
will enable us to see strong gravitational phenomena such as coalescence of compact stars 
or black holes. In order to extract signals from noise, we need to prepare in advance a list 
of expected templates of gravitational waves, which depend on many parameters of sources 
(masses and spins of compact objects, orbital angular momentum, eccentricity of their orbits, 
■ ■ ■). Therefore one of the current most important subjects in general relativity is to predict 
a precise waveform of gravitational waves for a set of given parameters. 

Furthermore, if we can directly observe gravitational waves, we may also put some con¬ 
straints on the theories of gravity, or we might End some evidence for alternative theories of 
gravity instead of general relativity. Although many alternatives have been proposed, most 
of them have been rejected by experiments or observations except for a few theories such 
as scalar tensor theories of gravity. As the simplest and proto-type among all scalar-tensor 
theories of gravity, the Brans-Dicke (BD) theory of gravity is often considered as an 
alternative. We know the strictest bound for BD theory is to the BD parameter ujbd ^ 500 
(the inhnity limit oi ubd is just general relativity), which was made by the observation of 
the Shapiro time delay using Viking on Mars [^. The possibilities of testing gravity theory 
using gravitational radiation as a tool were also considered in the early days §0. but it was 
not effective in the weak gravitational held of the solar-system. 

Using the binary pulsar PSR1913-I-16, the consistency of general relativity was tested, 
as we mentioned above. This binary was also considered as a tool for checking other gravity 
theories [§|,0. Recent analysis by Will and Zaglauer [^, however, concluded that this source 
cannot be used for testing BD theory effectively. This is because the two neutron stars in 
the binary are believed to have nearly equal masses, and the limit is also sensitive to the 
star models, hence the restriction on ujbd was weak. These classical tests, however, are 
carried out only in the weak gravitational held. On the contrary, observing gravitational 
waves directly from compact sources such as a BH-BH system will clarify some phenomena 
in strong gravitational helds and may determine which theory of gravity is correct. In 
order to analyze it from observational data, we also have to prepare a list of templates 
of waveforms for other theories of gravity. In BD theory, gravitational waves have three 
modes, i.e., a scalar mode, which we call a scalar gravitational wave (SGW), as well as 
two tensor modes (-1- and x modes) of conventional gravitational waves, which we shall call 
tensor gravitational waves (TGWs) here. Will |]ll| studied a constrained bound of ujbd from 


observing the inspiralling phase of compact binaries, and concluded that a bound strongly 
depends on mass diherence between components. To hnd the expected templates of SGWs, 
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Shibata, Nakao and Nakamura (SNN) [|T^ calculated the waveform of SGWs from spherically 
symmetrical dust fluid collapse and concluded that the advanced LIGO may detect SGWs 
even if the BD coupling constant ojbd is 10'^ ~ 10®. 

In this paper, we extend SNN’s studies and hgure out more aspects of both SGWs and 
TGWs behavior in BD theory. Since a spherical symmetrical system is just an idealized 
model, if we consider a more realistic situation such as the case that spherically symmetry 
is no longer valid, we may wonder whether SGWs can be really detected alongside TGW 
contributions. If we can observe the signals with the full polarization information of gravita¬ 
tional waves, then we can extract the scalar modes and do not need to worry about whether 
the TGWs will dominate the SGWs. However, if we have few observational bases, we may 
not be able to separate the polarizations. Hence it may be important to clarify in which 
situation the SGWs can dominate. Therefore, we study gravitational waves in BD theory 
and compare those scalar and tensor modes, analyzing test particle motions around a rotat¬ 
ing black hole. We analyze two cases: First, we calculate the waveforms and emitted energy 
both of SGWs and of TGWs by a test particle plunging with a constant azimuthal angle 
into a rotating black hole. This linearized analysis using a test particle is known to have a 
reasonably good agreement with fully relativistic numerical simulations in some cases such 
as a head-on collision of two black holes |^. Our case, we believe, also mimics a collision 
of two compact objects in BD theory. 

As for gravitational waves from gravitational collapse such as a supernova explosion 
followed by a formation of a neutron star, we know that no tensor modes are generated 
in a spherically symmetric case. Although the SGWs are still emitted in such a case and 
the theory of gravity might be determined by observing such SGWs, an ansatz of spherical 
symmetry is extremely idealized. Therefore, we next analyze a non-spherical dust shell 
collapsing into a rotating black hole. Gomparing the SGWs and TGWs, we can see how 
much deviation from spherical symmetry still gives us the domination of the SGWs. We 
then discuss the possibility of identihcation of the SGWs without observing polarization. 

This paper is organized as follows. In §2, we present our basic equations for calculating 
gravitational radiation in BD theory. We then analyze the waveforms and energy both of 
the tensor and of the scalar modes by a test particle plunging into a rotating black hole in 
§3. We also analyze non-spherically symmetric dust-shell collapse onto a Kerr black hole in 
§4. §5 is devoted to our conclusions and some remarks. Throughout this paper, we use the 
units of c = Gat = 1 |]I^ and the notations and dehnitions such as Ghristoffel symbols and 
curvature follow Misner-Thorne-Wheeler [113 . 


II. BRANS-DICKE THEORY AND BLACK HOLE PERTURBATIONS 

BD theory [Q is the simplest version of scalar-tensor theories of gravity in which a scalar 
held couples with metric tensor helds. The action S of BD theory is given by 


S= 


Ibvr 




+ s, 


matter 1 


( 1 ) 


where the coupling constant ujbd is called the BD parameter. The theory is reduced to 
general relativity in the limit of > cx). The conservative observational bound on ujbd 
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is u!bd > 500, which is obtained from the Shapiro time delay experiment 
From the action (|l|), the field equations become 


r - T -u 
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where = Rfiu — ^dfiuR is the Einstein tensor. Before we discuss perturbations of a black 
hole spacetime, it may be more convenient to introduce the Einstein frame, in which the 
metric and scalar perturbations are decoupled. Using the conformal transformation from 
the physical (Jordan-Brans-Dicke) frame to the Einstein frame as 
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fliy 
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the basic equations (j^) and are written as 
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where 0o is a constant, which is fixed by our ansatz Gat = 1 
the Einstein frame defined by 
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and $ is a scalar field in 
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From these equations, it is easy to see that any vacuum solutions in Einstein theory are 
equivalent with those in BD theory with <F= constant. In particular, the Kerr solution with 
<h = 0 is a solution in BD theory as well. 

We now discuss gravitational waves emitted by a test particle plunging into a Kerr black 
hole. In order to calculate the waveform of emitted gravitational radiation and its energy, 
we perturb the Kerr background spacetime. Since we are interested in a test particle motion, 
we adopt the matter Lagrangian as 
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where r and /i are proper time and mass of a test particle, respectively. In the case that 
the self-gravity of the test particle cannot be ignored, for example an extrapolation of a test 
particle into a nentron star or a black hole, /x shonld depend on the BD scalar 0, which 
does not gnarantee that the test particle moves along the geodesic and the gravitational 
weak eqnivalence principle is no longer valid (Nordtvert effect). In this paper, however, we 
consider only the case of /x=constant just for simplicity. When we extrapolate our results 
into the case of the direct collision of two compact objects, we should take into account the 
0-dependence in /i [|], but the difference may not be so large unless the test particle is a 
black hole, in which case scalar gravitational waves will not be emitted, at least within the 

T7[| . The energy-momentum tensor of the test particle is 


linear perturbation approximation 
now 


where z^{t) is the geodesic of a test particle. 

Now we linearize our basic equations in the Einstein frame around a background vacuum 
spacetime : 


glFJ = si”’ + h? “d 4 = 4'») + 4<‘), 

where = 0. Then the hrst order perturbation equations are 

GW[Ai?’l = = 0] 


where 
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is the hrst order perturbation of the Einstein tensor G^i, and denotes the D’Alembertian 
of the background spacetime. Since mass of the test particle /i is a hrst order quantity, both 
energy-momentum tensors in Einstein frame and in Jordan-Brans-Dicke frame diher only by 
a constant factor l/0o in the linear perturbation equations (|^), ([T5|) . We also hnd that the 
metric perturbation equation ( 0 ) for hlfj is exactly the same as those in Einstein theory 

0. 

In the Jordan-Brans-Dicke frame, the perturbation variables, which are expanded as 


9tiy = gf}+ and 0 = 0o + 

are given by the variables in the Einstein frame as 


( 17 ) 
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For the equation for from Eq. ( p!5[) or (^, we obtain 

Stt 
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(19) 


As for the metric perturbations in the Jordan-Brans-Dicke frame, we just have to study the 
perturbations and gravitational waves in Einstein theory, since the basic equations in the 
Einstein frame are exactly the same as those in the Einstein theory, which we have better 
knowledge in treating in ordinary procedures. 


III. GRAVITATIONAL WAVES FROM TEST PARTICLES FALLING INTO 

KERR BLACK HOLE 

A. Metric Perturbations 


We start considering a test particle motion in the Kerr background spacetime and analyze 
hrst the tensor modes of its gravitational waves (TGWs) emitted by a plunging test particle. 
In order to use these results in the analysis in §4, we take into account a test particle not 
only moving in the equatorial plane but also moving off the equatorial plane, which results 
have not been obtained so far. This means that we will also present new results in Einstein 
theory, i.e., the waveforms and energy flux in the case of test particles falling with a constant 
azimuthal angle into a Kerr black hole. 

Let us begin from a brief review of a black hole perturbation analysis in Einstein theory, 
which gives the same results for tensor modes in BD theory up to a factor l/0o- In order 
to calculate metric perturbations of a Kerr spacetime in the Einstein theory, we adopt the 
Sasaki-Nakamura formalism [|T^], in which the Regge-Wheeler equation is generalized by a 
transformation from the Teukolsky equation [^. The radial wave equation becomes 
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where the tortoise coordinate r* is dehned as 
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( 21 ) 


with A = — 2Mr + a?. The potential functions F{r) and U{r) are given as 

A 1 dy 
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Vt is the potential in the Teukolsky equation, which is given as 


VT{r) = 


+ 4i(r - M)K 
A 


+ 8iwr + A, 


( 22 ) 


where K = {r"^ + aP)w — ma. In Appendix A, we present the explicit form of the source term 
Simujir) for the general motion of a test particle. 

To describe the wave function Ximuir) of Eq. (1^) using the Green function method, we 
need two independent homogeneous solutions of the Sasaki-Nakamura equation, i.e.. 
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where k = w — ma/2r+ with r+ = M + Vbeing the radius of the event horizon, 
and the Wronskian 
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For observation, we only have to know the asymptotic behavior of the wave function Ximwir) 
at inhnity, which is 
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Then the waveform at inhnity is 
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where sSZi^) is a spin-weighted spheroidal function, which obeys 
Id' 
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with s = —2, —1, 0,1 or 2 and A being the eigenvalue. gSZi^) is normalized as 

rd9 sin 9USZ{0)f = l. 

JO 

From Eq. (^5]), the energy flux of gravitational waves at inhnity becomes 


dEAG^) 

dt 


327r 


= -E 

TT , 
L,m 


■ dE^rp dE^rp 

~dt 

roo 

/ duj 
J —00 Cq 


(27) 


Then, the total energy and the energy spectrum of the TGW are given as 
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respectively. 

Solving this equation, Kojima and Nakamura m analyzed the gravitational radiation 
induced by a test particle with mass n and angular momentum plunging on the equatorial 
plane or along the rotation axis into a Kerr black hole with mass M and angular momentum 
J = Ma. They calculated the waveform (only mode) and the emitted energy. They 
pointed out that the behaviors are similar to those of a Schwarzschild black hole but the 
amplitude becomes larger as the rotation parameter a gets larger. They showed that the 
emitted energy from the particle moving along the axis is smaller than that from the particle 
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moving on the equatorial plane. They also calculated the total energy, linear momentum 
and angular momentum extracted from the system by gravitational radiation and showed 
that all of them become larger as the rotation of black hole gets faster. 

For TGWs in BD theory, we have to multiply l/0o by the energy momentum tensor, 
and then the energy flux. But since 1 < 0o < 1-001 ioi ujbd > 500, we may ignore such a 
factor here. In order to compare the tensor modes with the scalar mode, we calculate the 
waveform and emitted energy of gravitational radiation induced by a test particle plunging 
with a constant azimuthal angle into a Kerr black hole. The particle motion with a constant 
azimuthal angle is a geodesic in Kerr spacetime, ii E = ^ and = 0, which we assume here 
(see Appendix). We have analyzed such a special geodesic because in the next section we 
will use it in the analysis of non-spherical dust shell collapse . We wish to emphasize that 
the results with the azimuthal angle 0 7 ^ 0 or 7 ^ 7 r /2 are new even in Einstein theory 


We have checked our numerical codes with the latest results, Mino, Shibata and Tanaka 


We show the emitted waveforms and emitted energy for several values of the Kerr param¬ 
eter a (a = 0, 0.5, 0.9, 0.99) and of the azimuthal angle 6 {6 = tt/O, tt/S, 7 r /2 (the equatorial 
plane)) in Figs. 1 and 2. We assume that the observer is in the position oi 9 = 7 r/ 2 , 0 = 0. 
One of the important results is that only in the case of an observer perpendicular to the 
direction of a direct collision, and for very small rotation of compact objects, do we hnd 
little gravitational radiation (Fig. 1 a = 0 case). This is consistent with a full numerical 
analysis of two black hole collision [|T3 . 


We also present the total energy of TGWs in terms of the azimuthal angle in Fig. 2. 
The emitted energy changes smoothly in terms of the azimuthal angle from the largest one 
for 0 = 7r/2 to the smallest one for 0 = 0, each values are found in [^| and W^ - 


B. Scalar Perturbations 


Next we consider the scalar mode 0 ^) of the perturbation. We expand 0^1 using spin- 
weighted spheroidal function as 
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then, from the held equation ([l^ , the radial wave function Zimujif) obeys the differential 
equation 
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and the trace of energy momentum tensor is expanded as 
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with E = cos^ 6. We find that 


TimLj{r) 


1 

(27r)3/2 


dr^r) 

dr 




(34) 


where we have to insert a test particle orbit described in terms of r, i.e., = 

{t{r),r,6{r),ip{r)). We again use the Green function method to derive the wave function 
Defining two independent homogeneous solutions as 
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We only need the asymptotic behavior of the wave function Zimuj{r) at infinity for observa¬ 
tion, which is 
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Then the waveform at infinity is derived as 
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The energy flux of the SGW observed at infinity is defined as 0 
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Then the total energy is now 
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In order to check our numerical code, for the case of a test particle falling into 
Schwarzschild background, we have compared our results with those by Shibata, Nakao 
and Nakamura [^], finding a good agreement. We show our results in Figs. 3 to 6 for 
a=0, 0.5, 0.9, and 0.99 and for 6 = 7r/6,7r/3 and 7r/2. Here we assume ujbd = 500. Since 
ijJBD ^ 1, the amplitude and energy of the SGWs are both approximately proportional to 
^/^BD (see Eqs. (^ID and (^)) and then we can rescale our results for any arbitrary value 
of ujbd- 

From the waveforms (Fig. 3), energy spectrum (Fig.4) and the total energy (Fig. 5), we 
find that the results depend only slightly on the values of a and 6. This 6 independence may 
be understood from the fact that the s wave (/ = 0) is dominant in the SGWs. As for the 
Kerr parameter, although the phase difference depends on a, we find that the waveforms 
themselves little depends on a. The energy spectrum has three main peaks (see Fig. 4) 
and the frequency of the largest peak is almost the same as that of the quasi-normal mode 
(QNM). These QNMs of scalar modes (/ = 0) are obtained by Simone-Will [Q and Anders- 
son for the Schwarzschild black hole case, using WKB approximation and phase-integral 
method, respectively. Here, using Leaver’s technique la of continued fractions, we derive 
the I = 0 modes of QNMs for a Kerr black hole and try to fit the tail part of the obtained 
waveform by them. In Fig. 6, we show both the calculated waveform (bold line) and the 
curve fitted by QNM (dotted line). The results show that the ringing tail of the SGW can 
also fit by the QNM. To draw a fitting curve, we need only the first QNM, i.e.. 


= ^eimK](i-r*+b) cos(Re[a;i](f - r* + b)). 


(41) 


with A = -0.0012, b = -33, Mui = 0.11045 - 0.10490z for a = 0, A = 0.0013,5 = 
-39,Mcui = 0.11251 - O.lOOOOi for a = 0.5M, A = -0.0014, b = -72, Mui = 0.11380 - 
0.09160i for a = 0.9M. However,the QNM fit to the ringing tail of the SGW is not as good 
as the fit to the TGWs. The situation does not change even if we include the higher QNMs. 
This may be because we cannot separate the burst part and the ringing tail part of SGWs, 
we cannot fit the tail part of SGW only by the QNMs. 

Because of the ratio of the maximum amplitude of GW to that of SGW is 200:1, a 
detection of the SGW itself is usually very difficult, unless we observe the polarizations and 
divide the tensor and the scalar modes. Only one exceptional case is that the observer’s line 
of sight is in the direction of a falling test particle into a Schwarzschild black hole. In this 
case, the SGW always dominates the tensor modes. 
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IV. GRAVITATIONAL WAVES FROM NON-SPHERICAL DUST SHELL 

COLLAPSE 


As we found in the previous section, when a test particle is falling into a Kerr black hole, 
the ratio of maximum amplitude of SGW to TGW is 1:200. This suggest that direct detection 
of SGWs from a collision of compact bodies is a hard task. However, the fact that SGWs 
are generated even in a spherically symmetric system, while TGWs require a deviation from 
spherical symmetry, gives us an expectation that the SGW may become dominant in some 
situations such as an almost spherical supernova explosion. Since a completely spherically 
symmetric system is unrealistic, we may ask how much deviation is allowed for the SGW 
to dominate TGW. In this section, we analyze the gravitational radiation from a collapsing 
non-spherical dust-shell, which consists of test particles (so its self gravity is ignored), into a 
Kerr black hole. In early studies by Nakamura and Sasaki and Shapiro and Wasserman 


2^, the TGWs for a deformed shell falling into a Schwarzschild black hole were investigated. 


Our study is an extension of their researches to a rotating black hole case and/or to BD 
theory. 

The gravitational radiation from a dust shell is derived from a linear combination of test 
particles falling into a black hole. If the dust shell is spherically symmetric, a phase can¬ 
cellation mechanism reduces the energy of gravitational waves to zero in the Schwarzschild 
black hole case. 

Let us begin by dehning the shape of a dust shell by a function r = rj(6*) {i = 1, 2); 


ri{6) = ro(l -I- 5 cos^ 6) (<5 > 0) (prolate), 

= ro(l -I- 5 sin^ 9) {6 > 0) (oblate). 


(42) 

(43) 


where 5 is a parameter which denotes a deviation from a spherically symmetric dust shell 
and tq is a typical radius of the shell. 

In the case of a Schwarzschild black hole, test particle motion with zero angular mo¬ 
mentum does not depend on the azimuthal angle 6, and the dynamics of particles in the 
shell can be replaced by that on the equatorial plane. This fact makes our calculation of 
energy of gravitational waves easy as follows. Using the shape function ri{9), we introduce 
a correction function fi^Lj as 

fim. = l TdO Sin 9 (44) 

2 JO 


where Pi is the Legendre function and t{r) is the coordinate time when the particle comes to 
a position r. Then we dehne the source term Simco and wave function Ximuj for non-spherical 
dust collapse from those of test particles as 


[shell] n [particle] 

[shell] p -^[particle] 


(45) 

(46) 


We then obtain a total energy of gravitational radiation from non-spherical dust shell collapse 
as 
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i \ dui ) 


( 47 ) 


Imuj 


where 


/ rf^(TGW) \ 

V Jlmu. 


I flmuj I 


8w \flmu! 


V du 

X {oo)[particle] 
u 
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Co 


As for the SGW, we have 


fe; •'-“ \ doj j 


( 48 ) 


(49) 


where 


/^^(SGW)\ 
ha; j 


Imuj 


{2uJbD + 
16(a;s£) + 2 ) 


Imuj 


\r \2\/^{oo)[particle]\2 
\ JImuj \ \ ^Imu) \ * 


( 50 ) 


In a Kerr spacetime, however, we cannot adopt this simple formula because the radial 
wave function depends on the azimuthal angle 9, which means that we cannot factorize 


a correction function as Eq. (|48|). In this case, we have to calculate the total energy by 
summing up test particles. Therefore, as for test particles of the dust shell, we choose 
E = fi and Lz = 0, for which a trajectory of 6*=constant becomes a geodesic. First, using 
Eq. (^9|) , we prepare the wave function from a test particle with constant 6*, which we denote 
Ximu}{,'r] 6). Then we calculate the energy spectrum of the TGW as 


'rf^(TGW)\ 


duj 


■J 


= 2u^ 


Imuj 


fTT 

/ dOsmO^^ 

'0 Co 


r;9) 


^iw[t{ro)-t{ri{d))] 


•^mO? 


(51) 


where 9) is the asymptotic form of Ximcj{r] 9) at inhnity. To evaluate the integration 

(pl), we prepare data at 60 points for 9 in [0, ttI. The total energy is given by Eq. 

As for the SGW, we also End 

{2ujbd + 3 )^ 



Imuj 


16{u!bd + 2 ) 


XU 


d9sin99) e 


Jw[t{ro)-t{ri{9))] 


■'mO- 


(52) 


We show our results in Figs. 7 and 8. For a Schwarzschild background, Nakamura and 
Sasaki calculated the total energy of TGWs from a dust shell and pointed out that the 
energy will be generated maximally for a deviation of ~ 0.5 from spherical symmetry. 

The reason is that if the dust shell is close to a spherical symmetry, the TGW will not be 
generated so much while if the deviation is larger than a phase cancellation effect 

will also reduce the emitted energy . 
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For the case of SGW, we find that the energy does not depend on the deformation 
parameter <5 if <5 is smaller than some critical value = 1.5) (see Fig. 7). This is 

understood by the fact that the /=0 mode is dominant for the SGW. Beyond the critical 
value, the above cancellation effect also reduces the energy of the SGWs. At the critical 
deviation, the time lag between the longest and shortest radii is almost the same as the most 
dominant QNM period. If the deviation is larger than the critical value, i.e., if the time 
lag is larger than the most dominant QNM period, the total energy is reduced by a phase 
cancellation. In fact, the difference between the values of and is explained by 

the difference between QNM periods of TGWs and SGWs. 

For the case of a rotating black hole, we show the results in Fig. 8 . The results are 
similar to those of the Schwarzschild case except for the following two points. First, since 
the system is not completely spherically symmetric, the TGW is also generated for any values 
of S. There exists some value of (5o, which is not zero and depends on the rotation parameter 
a (e.g. (5o ~ 0.02 for a = 0.5, (5o ~ 0.08 for a = 0.9, and 5o ~ 0.10 for a = 0.99), at which 
the emitted energy of TGWs takes a minimum value. As a increases, its minimum energy 
and the corresponding deformation parameter (5o increase (see Fig. 8 (b)). The increase of 
minimum energy is easily understood because of a large deviation from spherical symmetry. 
The fact that (5o increases as a increases may also be explained by the following reason. 
When we estimate the falling time of a test particle from a hnite radius in a Kerr spacetime, 
the time along the rotation axis is longer than that on the equatorial plane if the initial radii 
are equal. Then if the dust shell is slightly oblate, the time lag between test particles along 
the rotation axis and on the equatorial plane becomes smaller, which is close to a spherically 
symmetric case, and then a phase cancellation effect will reduce the emitted gravitational 
radiation. As for the SGW in a Kerr geometry, the results are almost same as those in a 
Schwarzschild one. 

Secondly, the existence of a non-vanishing minimum value of TGW is important for 
observation. Since the SGW does not depend on 6 , we can hnd how much deviation is 
allowed for SGW to dominate TGW from Figs. 7 and 8. That is, the energy of the SGW is 
larger than that of the TGW if (5i < 5 < ^2 and ubd = 500 ((5i = —0.07, S 2 = 0.07 for a = 0, 
= -0.05, 62 = 0.09 for a = 0.5, (5i = 0.03, 82 = 0.11 for a = 0.9, and = 0.05, 62 = 0.12 
for a = 0.99 when ujbd = 500). The emitted SGW is proportional to 1/ojbd if ^bd 3> 1. 
Then if ojbd is larger than some critical value Ubbi which depends on a, the SGWs never 
dominate TGWs {uJbb = 20, 000 for a = 0.5, Ubb = 1, 250 for a = 0.9, and Ubb = 750 
for a = 0.99). We may speculate that in generic gravitational collapse the SGW does not 
dominate the TGWs if ubd > several thousand. 


V. CONCLUDING REMARKS 

In this paper, analyzing test particles around a Kerr black hole, we studied both scalar 
and tensor modes of gravitational waves in BD theory. We examined two models: A test 
particle plunging with a constant azimuthal angle into a Kerr black hole, and a non-spherical 
dust shell collapsing into a Kerr black hole. In the hrst model, we show that the SGWs have 
little dependence on the the rotation parameter a and the azimuthal angle 6 because the 
s-wave is dominant. If a coupling constant of BD theory ojbd ^ 500, however, the TGWs 
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dominate the SGWs except for one very limited case, that is, an observation in the face-on 
direction of a test particle falling into a Schwarzschild black hole. In the second model, 
we also show that the SGWs dominate only for the highly spherically symmetric case. In 
particular, for the case of a rotating black hole, the emitted energy for TGWs takes a non¬ 
zero minimum value when the deformation is slightly oblate. Since that for SGW is less 
sensitive to a deviation from a spherical symmetry, the SGW cannot be dominant unless 
^BD ^ 750 for a/M = 0.99 or unless ubd ^ 20, 000 for a/M = 0.5. 

Here, we remark on the detectability of SGWs. In the First/Advanced LIGO operation, 
the observation limit of the gravitational wave amplitude is expected to be h > 10“^^ for 
the 50 ~ 200 Hz band, and h > 10“^^ for the 10 ~ 1, 000 Hz band, respectively [Q. 

First we mention that the QNM frequency provides some naive measure for the observable 
range. The relationship between the observed frequency / and the QNM frequency uj is 

/ ~ 3.30 X 10^ Mcu [Hz]. (53) 

for typical observation of NS (M ~ I. 4 M 0 ) and BH (M ~ IOMq). The most dominant 
contribution of QNMs is n = 1, for which we hnd Re Moo ~ 0.37 for TGW and Re Moo ~ 0.11 
for SGW. Hence, for the BH-NS system, the observable frequencies in First (Advanced) 
LIGO are Moo = 0.02 ~ 0.08 (0.004 ~ 0.4), while for the BH-BH system, 0.15 ~ 0.61(0.03 ~ 
3.04). Both observations of SGW and TGW are within preferable ranges of frequency. 

Secondly, we estimate the maximum amplitude of gravitational radiation. The expected 
amplitudes of SGW and TGWs (h^^GW)^ are 


fc ( TGW , ^ 2 ,, (^) ( HEH ) . 


where we set = 0.4 and = 0.002 from our results. 

Fixing the BD coupling constant oobd = 500, we hnd that the limiting distance of SGWs 
in First LIGO for NS-NS, BH-BH collision is r=100kpc, while r=20Mpc for TGWs. This 
limiting distance becomes much longer when we have the NS-BH system and Advanced 
LIGO detector, i.e., the target for SGW could set farther than the Virgo Gluster (20Mpc). 

Gonversely, if the SGW is observed in a gravitational wave detector with polarization 
information, what kind of constraint on the BD coupling constant we will hnd ? We 
can at least discuss the upper bound. For example, we can easily say that if NS-BH collision 
occurred in our Galaxy and the SGW were observed, we could determine dehnitely the BD 
coupling constant even for much higher values than 500 (see Table 1). 

Although our analysis has been done by test particles, we believe that our conclusions 
for the SGWs are valid for more generic cases. As for other scalar-tensor theories of gravity, 
we expect similar results, although we may need further analysis order to determine 


a theory of gravity by gravitational wave detection, the above results require observations 
of polarizations and extractions of information of the scalar mode from observed data, for 
which we need several observational sites in the world. 

It is very interesting to test theories of gravity by direct observation of gravitational 
radiation. It will open a new window in gravitational astronomy, which may provide a third 
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and the final observational eye to see the Universe and may be expected to reveal some 
problems in fundamental physics (such as the equation of state at high density) as well as 
to determine many properties and orbital elements of compact objects. 
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APPENDIX A: SOURCE TERM OF GRAVITATIONAL WAVES BY A TEST 

PARTICLE IN A KERR SPACETIME 


In §. 2, we summarize the tensor gravitational waves from a test particle in a Kerr 
spacetime. In order to solve the basic equation (the Sasaki-Nakamura equation), we hrst 
have to give the source term Simuj in Eq. (pOl). Nakamura and Sasaki gave the explicit form 
for a test particle moving on the rotation axis |l^ and Kojima and Nakamura gave it on the 
equatorial plane |^. Shibata presented it for a motion with r = constant [^. However, 


no one has so far written it down for generic motions. Since we had to calculate the case 
with off-plane motion {6 = constant), here we present the explicit form of a source term for 
generic motions in use of numerical calculation. Here we consider only non-periodic motion 
or unbound system. If the system is bounded, we may hnd many turning points, for which 
our expression is no longer valid. We have to rewrite the orbital motion in terms of time 
(either a coordinate time f or a proper time r) instead of r. This may be straightforward 


A test particle in a Kerr spacetime is described by the equations of motion as 


where 


= ±VR, 

dr 

(Al) 

= ±v^, 
dr 

(A2) 


(A3) 

df 

S— = —a{aE sin^ 9 — L^) H-r—P, 

dr A 

(A4) 

P = E{r‘^ + a^) — aLz, 

(AS) 

P = - A (r^ + {Lz - aEf + c) , 

(A6) 

0 = C — cos^ 9 — F^) H——, 

V sin 9 1 

(A7) 


and E and are conserved energy and z-component angular momentum of a test particle, 
respectively, and C is the Carter constant. We can easily see that 6 = constant is a geodesic 
if E = /i and = 0. 

The source term Simui is given by 

yA U K 

J ^dr). (A8) 

Here W is divided into three parts as 


W — yVnn + y^mn + EV’mm; (^9) 

roo roo roc 

>V„n =/o(r)e*^Pi-1- / dr7i(r')e'^P9 + / dr" (AlO) 

Jr Jr Jr' 

y^fhn = ^o(r)e*>^Pi -F r dr'gi{r')jy^'^ -f H d^H dr"g 2 {r")jy^''\ (All) 

Jr Jr Jr' 

poc poo poo 

= ho(r)e*^Pi + / dr'h{r')J^y'^ + dr' dr%(r'0e*^P"i, (A12) 

Jr Jr Jr' 
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where 


X = uj(t + r*) — m(p, 

fo 9^nn) 

2 (dw. 


dip 


fi =-—+ i{ 2 ausm 6 - m)—WnnU 


f 2 = —. 




d Wnn ^ 


dr"^ 
+ 2 iaoj sin^ 6 


2 im 


dr dr 


2 2 ■ 4 n I \ 

-a u sin 0 \ — \ m 


dWnn dip d^ip 

dr dr dr"^ 

2 



Wr, 


dr 


_ ^ ( 1 ) ( 2 ) 
90 P '^mn 

OJ 


I ( 1 ) 
9l = -U!kn 
u 


_y;( 3 ) + y;( 2 ) ^ + 2 o — 

^mn ' ^mn 7 ' 7 

dr dr 


( 2 ) 


_ ^ ( 1 ) 
92 'UJmn 

OJ 


A. 

dr 


(3) dwk 


Wkn + P- 


dr 


-i (m — aoj sin^ 9^ wj^l — p 


(3) dwA^n \ dip 


dr / dr 


2^2 


ho — — 

hi = - 


1 r p- 


2 \dr/dT 

1 PP^ 

2 \dr/dT\ 


qauj („.,(}) \2 

—2^lm\^mn) j 





3 2' 

P r , 


with 


h - ^ PA C«‘-’('.,,(h'|2_^ / -4_^/ 3 2'|\ 

/ r Q 

-^dr, 


CLWmn ■ f • 2 n\ 

+ Ip —^- %[m — aujsm d] p w. 

dr V ' 

( 2 ) 


( 2 ) ^ 
fhn~i 

dr 


p = 


'^nn — 


r — ia cos 9 ’ 


pr 


2pYk 


dr- 


dr 


c\[f.-*c\uszpX 


^ V© + ismO{aE - 

Sin 9 

(_ 2 ) _ ^ P r\ ( qauj\ 

^mn 2 ^2{PP—2^lm) ■) 


(A13) 

(A14) 

(A15) 


(A16) 

(A17) 

(A18) 


(A19) 

(A 20 ) 

(A 21 ) 

(A 22 ) 

(A23) 

(A24) 

(A25) 

(A26) 

(A27) 


17 



[p^-■2SZ^{P^P , 

where an overbar denotes the complex conjugate. The operator is dehned by 

X (9 m 

= tt;: + ao; sin 6*-- + s cot 0 . 

oO sin 0 


(A28) 


(A29) 
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Figure Captions 


Fig. 1: Waveforms of TGWs by a test particle with constant azimuthal 
angle 6 falling into a Kerr black hole. We show two modes (/i+ and 
/lx) of TGWs for 6 = 7r/6 ((a), (b)), = 7r/3 ((c), (d)) and mode for 
9 = t:I2 ((e)) (hx vanishes.). We choose a Kerr parameter a= 0 (solid 
line), 0.5 (dashed line), 0.9 (dash-dotted line) and 0.99 (dotted line). 


Fig. 2: Emitted energy of TGWs showen in Fig. 1. 

Gircle (filled), square, triangle and circle (open) represents the cases of 
a = 0, a = 0.5, a = 0.9 and a = 0.99, respectively. 

Fig. 3; Waveforms of SGWs by a test particle with constant azimuthal 
angle 9 falling into a Kerr black hole for ujbd = 500. We show SGWs for 
0 = 7r/6 ((a)), = 7r/3 ((b)) and = 7r/2 ((c)). We choose a Kerr parameter 
a= 0 (solid line), 0.5 (dashed line) and 0.9 (dash-dotted line). We find 
that they do not depend very much on the azimuthal angle and the Kerr 
parameter, although we find a phase difference depends on a. 


Fig. 4: Energy spectrum of SGWs shown in Fig. 3. 

Solid, dashed, dash-dotted and dotted lines represent a = 0, a = 0.5, 
a = 0.9 and a = 0.99, respectively. We find that the frequencies which 
emittes largest energy spectrum corresponds to that of the QNM. All 
lines (solid line (a = 0), dashed one (a = 0.5), dash-dotted one (a = 0.9) 
and dotted one (a = 0.99)) almost coincide with each other. 


Fig. 5: Emitted energy of SGWs shown in Fig. 3. 

Gircle (hlled) , square, triangle and circle (open) represent a = 0, a = 
0.5, a = 0.9 and a = 0.99, respectively. 

They also do not depend very much on the azimuthal angle. 


Fig. 6; The tail part of the SGW is fitted by one of the QNMs, which 
corresponds to the largest peak in the energy spectrum (Fig. 4). Solid 
line represents the SGW, while dashed line denotes the waveform fitted 
by one of the QNMs ((a) a = 0, (b) a = 0.5 and (c) a = 0.9). 
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Fig. 7: Emitted energy by a non-spherical dust shell collapsing into 
a Schwarzschild black hole ((a) prolate shell and (b) oblate shell) for 
= 500. A spherically symmetric shell does not generate TGWs, 
while SGWs are emitted. Since the I = 0 mode dominates in the SGW, 
the emitted energy of the SGW does not depend on a deformation pa¬ 
rameter S. The SGW dominates the TGWs for (5i(= —0.07) < S < 52{= 
0.07). 


Fig. 8: Emitted energy by non-spherical dust shell collapsing into a 
Kerr black hole ((a) prolate shell and (b) oblate shell) foi ivbd = 500. 
Dashed, dash-dotted and dotted lines represent a = 0.5, a = 0.9 and a = 
0.99 respectively. Since the system is no longer spherically symmetric for 
any values of 5, we expect the TGWs as well as the SGW. At 5o, which 
is not zero and depends on a rotation parameter a (Jq ~ 0 .02, 0.08 and 
0.10 for a = 0.5, 0.9 and 0.99, respectively), the emitted energy of TGW 
takes a minimum value. As a increases, its minimum energy and the 
corresponding deformation parameter Sq increase. Since the I = 0 mode 
dominates in the SGW, the emitted energy of SGW does not depend on 
a deformation parameter 5. It also depends little on a. Both energies 
of TGW and SGW decrease for 6 ^ (5cr and 5 ^ because of 

a phase cancellation effect. The SGW dominates the TGWs for 5i(= 
-0.05, 0.03 and 0.05 for a = 0.5, 0.9 and 0.99) < (5 < 52 i= 0.09, 0.11 
and 0.12 for a = 0.5, 0.9 and 0.99). 
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TABLES 


Binary System 

Coma cluster 
(~ SOMpc) 

Virgo cluster 
(~ 20Mpc) 

Our Galaxy 
(~ lOkpc) 

NS-NS, BH-BH 

^BD ^ 60 

^ 250 

^ 5 xlO® 

NS-BH 

^ 400 

^ 1,700 

^ 3.5 xlO® 


Table 1: The upper bound of the BD coupling constant ujbd if the SGW with a polar¬ 

ization information independently of the TGW is detected. 
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